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A HIGHER BOLTZMANN DISTRIBUTION 


MICHAEL J. CATANZARO, VLADIMIR Y. CHERNYAK, 
AND JOHN R. KLEIN 


Abstract. We characterize the classical Boltzmann distribution 
as the unique solution to a certain combinatorial Hodge theory 
problem in homological degree zero on a finite graph. By sub¬ 
stituting for the graph a CW complex X and a choice of degree 
d < dimX, we define by direct analogy a higher dimensional Boltz¬ 
mann distribution as a certain real-valned cellular (d— l)-cycle. 
We then give an explicit formula for . 

We explain how these ideas relate to the Higher Kirchhoff Net¬ 
work Theorem of [CCK]. We also give an improved version of the 
Higher Matrix-Tree Theorems of [CCK]. 
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1. Introduction 

1.1. Motivation. Physics and chemistry are rife with processes in 
which some quantity varies with time in a complicated, irregular way. 
The evolution of such a system is often modeled by a so-called master 
equation, which takes the form p = Hp. Here, p is a time dependent 
distribution on the states of the system and H, the master operator. 
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governs time evolution. The stationary distribution of a master equa¬ 
tion is known as the Boltzmann distribution (or Gibbs measure). 

For example, the classical Boltzmann distribution from thermody¬ 
namics governs how the system’s states are populated, depending on 
their energy. Specifically, the probability of the system to be in state 
j with energy Ej is proportional to where jd = T is the 

temperature and fcs is the Boltzmann constant. The Boltzmann dis¬ 
tribution is then the normalized distribution 


P 


B 


Ej 


A Markov process with discrete state space is conveniently described 
by a state diagram. The latter is a directed graph T whose vertices label 
the states and whose edges label the transition probabilities between 
states. Each directed edge e of T is equipped with a transition rate k^. 

A certain class of processes is obtained by the following procedure: 
start with a connected undirected graph X (i.e., a CW complex of 
dimension one) with vertex set Xq and edge set Xi. Then T will be the 
double of X, i.e., the directed graph having the same vertices, where 
each edge of X is now replaced by a pair of opposing directed edges. 
Each directed edge of T is given by a pair {i, a), where Ms a vertex of 
X and a is an edge of X that is incident to i. The transition rates of 
the process are given as follows: choose a real number fd > 0 (inverse 
temperature) and functions E: Xq —)■ R and W: Xi —)■ R. Then the 
transition rate across {i, a) is given by 


^ 


When the transition rates are written in this way, the process is said 
be in Arrhenius form. 

The process just described is not completely general since T is a 
double. In fact, the above is an example of a process that is in detailed 
balance (or time reversible) in the sense that there exists a distribution 
TT: To —)■ R+ such that 


( 1 ) E{i)ki^oj = 

for any edge a of X with joining a pair of vertices i and j (see e.g., [Ke, 
ch. 1]; in our example, take 7r(i) = The number nii^ki^a is the 

probability flux across (qa). Equation (1) says that the pair of states 
i and j are in equilibrium along a. It is not difficult to show that every 
process (with discrete time and hnite state space) in detailed balance 
can be written in Arrhenius form. 
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If X is finite, the Boltzmann distribution arises as the normalized 
vector 

Z:='^7r(i). 
ieXo ieXo 

Consequently, the Boltzmann distribution is a normalized equilibrium 
distribution for the process. 

The higher Boltzmann distribution developed in this paper is asso¬ 
ciated with a Markov process on a CW complex which is a natural 
generalization of the process we described above on a graph. With 
the intention of providing motivation, the process is described below. 
However, we will not investigate its properties here as is not within the 
current scope. 

Let X be a hnite connected CW complex, let d < dimX be a hxed 
positive integer, and let denote the set of /c-cells. With respect 
to a mild technical assumption, we will show how associate a Markov 
process. A state of the process will be an integer-valued cellular (d —1)- 
cycle on X in a hxed homology class. Roughly, a directed edge from 
( to another state (' is given by an “elementary homology” in the sense 
that 

C = C + ude, 

where u is an integer and e is a d-cell that is incident to ^ at a specihed 
(d — l)-cell /. We also require that e not be incident to (' at /. 

When d = dimX = 1, it will turn out that the process coincides 
with the one constructed above. Before giving the details, we explain 
the technical assumption on X, which amounts a homological condition 
that will eventually be removed as it is not required for the results of 
this paper: 

Definition 1.1. The CW complex X is d-pseudo-regular ii for every d 
cell e G Xd and any (d — l)-cell / G X^-i, we have 

{dej) G {0,±1}, 

where the above denotes the coefficient of / in the boundary of e. If 
X is regular then it is d-pseudo-regular for all d. In particular any 
connected polytopal complex [Z, chap. 5] or any connected hnite sim- 
plicial complex is d-pseudo-regular. If d = dimX = 1, then X is 
automatically d-pseudo-regular. 

Assuming X is d-pseudo-regular, we can now describe the process. 
Consider the oriented graph in which a vertex is given by an integer 
(d — l)-cycle C, G Zd_i(X;Z). A directed edge from a vertex to a 
vertex C,' is given by a pair 
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in which 

• e G Xd is a (i-cell and / G Xd-i is a (d — l)-cell; 

• (C,/) 7^ 0 7^ {dej) ; 

• (C',/) = 0. 

Note that the above conditions imply 

C' = C - {CJ){deJ)de, 

so C is obtained from ( by means of an elementary homology. We 
denote the directed edge by (/, e, (). 

We now fix a base vertex 

X G Zd-i{X; Z) , 

such that the homology class [f] G Hd-i{X] Z) is non-triviah Consider 
the subgraph generated by x in the sense that a vertex ( lies in this 
subgraph if there exists a finite sequence of directed edges x = (o ^ 
Cl Cfc = C) i-G-, a finite directed path from x to (. An edge 

belongs to this subgraph if it occurs in such a path. We denote this 
subgraph by 

^X,x ■ 

Definition 1.2. The cycle incidence graph of {X,x) is the directed 
graph Tx,x. 

Example 1.3. Suppose X has dimension one and a; := i is any vertex 
of X. In this situation T^,* coincides with the double of X. 

Remark 1.4. If dimX > 1 then the cycle-incidence graph is usually not 
a double, so the process won’t be detailed balance. Furthermore, Fx,* 
typically has infinitely many vertices. 

Example 1.5. Let X be the two dimensional torus with regular cell 
structure indicated by the following picture: 


ei 

62 

es 

64 


Then X is obtained from the picture by identifying the opposite sides 
of the outer square. Consequently, there are 4 two-cells (ci,..., 64 ), 8 
one-cells and 4 zero-cells. Our initial one-cycle Co = ^ is fhe meridian 
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given by the red segment. In this case the cycle incidence graph has 
inhnitely many vertices. To see this start at x and jump across Ci. 
Modulo a choice of orientations for the cells, the jump results in the 
cycle given by Ci = x+dei, and the latter is incident to the 62 . Jumping 
across 62 results in the cycle C 2 = x + de\ — which is incident 
to Cl. Moreover C 2 is distinct from x. Iterating this procedure by 
jumping hrst across e\ and then across 62 , results in an inhnite number 
of distinct vertices ( 0 , ( 1 ,( 2 , ■■ ■ (with ( 2 k+i = x + (2k + l)dei — ( 2 k)de 2 
and ( 2 k = X + (2k)(dei — de 2 )) in the cycle incidence graph of (X,x). 


To complete the description of the process, we must label the edges 
of Tx,x with transition rates. This is done as follows: choose functions 
E: Xd-i —)■ M and W: X^ —)■ M and an inverse temperature /9 > 0. 
Then the transition rate of the directed edge (/, e, () is given by 

The description of the process is now complete. 

Summarizing, as the initial (d — l)-cycle evolves with respect to the 
process, it jumps from cycle to cycle by adding/subtracting the bound¬ 
aries of d-cells, analogous to a particle ‘jumping’ over an edge on a 
graph. That is, the evolving cycle does not change its homology class, 
and so each homology class will have its own distinct dynamics. 

One can interpret the higher Boltzmann distribution dehned in this 
paper as the long-time limit of the average distribution within each 
homology class. By “average,” we mean the average over the stochastic 
process in the sense of probability theory. This will be pinned-down in 
a subsequent paper. 


1.2. The combinatorial Hodge problem. We now turn to the prob¬ 
lem of formulating the higher Boltzmann distribution in terms of com¬ 
binatorial Hodge theory. Let X be a hnite CW complex and let 
d < dimX be as above. Henceforth, we do not require X to be d- 
pseudo-regular. Let be the cellular chain complex of X with 

real coefficients. Each Cj{X-, M) is equipped with a standard inner prod¬ 
uct (—, —) which is determined by declaring Xj to be an orthonormal 
basis. 


Definition 1.6. A scalar structure on X consists of functions 
Ej Xj —y M, j = 0 , 1 ,... 

When j = d we typically write E := Ed-i and W := Ed. 
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Let /3 > 0 be given. Then the function Ej equips with a 

modified inner product 

{x, y)Ej := y) x,y e Xj . 

In our next formulation, we only make use of the function E = E^-i- 
Dehne the formal adjoint 

d*E: Cd-i{X-R) ^ CfiX-R) 

using the standard inner product on C'rf(X;R) and the modihed inner 
product on Cd-i(X; R), i.e., d’^ = d*e^^, where d* is the formal adjoint 
in the standard inner product structures. 

Definition 1.7. Let X and E be as above. The combinatorial Hodge 
problem in degree d — 1 is the following: given x G Hd-i{X]R), hnd an 
explicit formula for the unique cycle p G Zc;_i(X;R) such that 

• p represents x, and 

• p is co-closed, i.e., d*Ep = 0. 

The condition that p be co-closed can be re-stated as the assertion 
that p should be orthogonal to any boundary with respect to the mod¬ 
ihed inner product: 

{a,d%p) = {da,p)E = 0, 

for any a G X^. 

Remark 1.8. The combinatorial Hodge problem for (X, E) is equivalent 
to hnding an orthogonal splitting of the quotient homomorphism p 
appearing in the short exact sequence 

0 5,_i(X;R) Z,_i(X;R) ^ H,_fiX-R) 0, 

with respect to the modihed inner product on Zc;_i(X; R) C Crf_i(X; R). 

The original Hodge problem asks to hnd a unique harmonic represen¬ 
tative for any cohomology class on a compact, orientable Riemannian 
manifold. By relaxing the hypotheses to a connected CW complex, we 
are able to write down an explicit formula. Our solution to the com¬ 
binatorial Hodge problem will involve a summation over spanning co¬ 
trees. The latter are certain subcomplexes of X of dimension d—1 which 
are a higher dimensional analog of the vertices of a graph. They are ho¬ 
mologically dual to the higher dimensional spanning trees of [CCK] and 
hence their name. In the following dehnition, fik{X) = dimiLfc(X;Q) 
denotes the kth Betti number of X, and X^^^ denotes the fc-skeleton of 
X. 
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Definition 1.9. A {d — l)-dimensional spanning co-tree for X is a 
sub complex L G X such that 

(a) The inclusion L G X induces an isomorphism 

ie* ■ Q) ^ Hd-i{X; Q ); 

(b) (3d-2{L) = (3d-2{Xy, 

(c) gLg 

Clearly, the number of spanning co-trees is hnite. 

Remark 1.10. Equivalently, conditions (a)-(c) are equivalent to condi¬ 
tions (a’),(b) and (c), where 

(a’) The relative homology group Hd-i{X, L]Q) is trivial. 

Remark 1.11. We are indebted to a referee for pointing out that our 
notion of spanning co-tree is equivalent to that of a relatively acyclic 
complex in [DKM2] as well as to the complement of a cobase in [L, 
p. 156]. 

Spanning co-trees come packaged with auxiliary data that will be 
used for obtaining the desired splitting. Observe that the projection 
Zd-i{L] Z) —)■ Hd-i{L] Z) is an isomorphism since L has no d-cells. Let 
0L be the composite 

(3) (Pl: Zd-i{L-Z) ^ Hd-i{L-Z) ^ Hd-i{X-Z). 

Then (pi becomes an isomorphism after tensoring with the rational 
numbers by the dehning properties of L. We invert (pi rationally to 
obtain a homomorphism of rational vector spaces 

(4) 'iPi-. Hd-i{X- Q) Zd-i{L- Q) ^ Zd-i{X- Q). 

Since Hd-i{X, L-,<Q) is trivial, the group Hd-i{X,L]7j) is finite. Let 

(5) aL-.= \Hd-i{X,L-,Z)\ 
denote its order. 

We define the weight of L to be the real number 

(6) n = tl{E) :=al JJ . 

With respect to the above definitions, we can state the solution to 
the combinatorial Hodge problem: 
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Theorem A (Boltzmann Splitting Formula). Let {X,E) and x G 
be as above. Then the solution to the combinatorial Hodge 
problem is given by p = in which : Hd-i{X-,M.) —)■ Zrf_i(X;M) 

is the homomorphism 

L 

where the sum runs over all {d — 1)-dimensional spanning co-trees L 
and V = 

Theorem A enables us to define the higher Boltzmann distribution: 

Definition 1.12. Let {X,E) be as above and let x G Hd-i{X-,7j) be 
an integer homology class. The higher Boltzmann distribution at x is 
the real {d — l)-cycle 

L 

where x G Hd-i{X] Q) is the image of x with respect to the homomor¬ 
phism Hd-i{X]Z) -)■ Hd-i{X]Q). 

Remark 1.13. The classical Boltzmann distribution is a special case 
of Definition 1.12: let X be finite connected graph and take d = 1. 
Then the spanning co-trees of X are given by the vertices. We take 
X G Hq{X] Z) = Z to be the canonical generator (given by choosing a 
vertex of X ; the generator is independent of this choice). For a vertex 
L = j the normalized weight is given by 

VVi = , Z =Y , 

ieXo 

since (pL is sji integral isomorphism. Then 'ipiix) = j and the assertion 
follows. 

Remark 1.14. When = 0, the coefficients tl are rational numbers 
and the map T is defined as a homomorphism of rational vector spaces 
Hd-i{X]Q) —)■ Zd-i{X-,Q). If we further assume that x is a rational 
homology class, then the solution to the combinatorial Hodge problem 
gives an explicit expression for the Harmonic “forms” with respect to 
the combinatorial Laplacian —{dd*-\-d*d): C'rf_i(X; Q) —)■ (7^-1 (X; Q). 

Remark 1.15. The proof we give of Theorem A is an application of the 
theory of generalized inverses to the quotient map p\ Zd_i(X;M) —)■ 
Ffrf_i(X;M) (cf. [M], [P], [BG]). In the late 1980s a summation formula 
was given for the Moore-Penrose pseudo-inverse [Be], [BTj; this is the 
formula we make use of. When writing the current paper, we also came 
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to realize that by applying the summation formula to split the inclusion 
map 

one gets another proof of our higher dimensional analog of Kirchhoff’s 
theorem on electrical networks [CCK] (see Remark 3.6 below). 

Remark 1.16. The main application of Theorem A will appear in the 
hrst author’s Ph. D. thesis [C]. 

Remark 1.17. It is tempting to speculate whether a result like Theo¬ 
rem A holds in the case of a Riemannian manifold. For this, one needs 
notion of spanning co-tree adapted to the space of differential forms. 
Then the sum appearing in the Theorem A would presumably be re¬ 
placed by a convergent inhnite series of operators indexed over the set 
of spanning co-trees. This would give an explicit solution to the main 
result of classical Hodge theory. 


1.3. Asymptotic behavior. Given {X,E) and d as above, with E 
suitably generic, it turns out that sum of Theorem A is asymptotic as 
a function of [3 to the term with highest weight. To explain this, let 
E* := denote the set of (d — l)-dimensional spanning co-trees 

ofX. Let 

E^, ■. E* 

be the functional given by 

L I—)■ Eh ■ 

b&Ld-x 


A scalar function E\ Xd-i —)■ M is said to be non-degenerate if it is 
one-to-one. This is clearly a generic condition. 

As spanning co-trees are a matroid basis [L, p. 156], a greedy algo¬ 
rithm shows that if E is non-degenerate then the function Ejr* possesses 
a unique minimum L^} 

As the parameter [3 tends to cxo, it is easy to check that the operator 
T appearing in Theorem A is asymptotic to when we consider 
these as vector-valued functions with components indexed over E* (see 
the proof of [CKS, lem. 3.6]). More precisely, we have 

Corollary B. Let {X,E) and x G Hd-i{X;h) be as above. Assume in 
addition that E : X^-i —)■ M zs non-degenerate. Then we have 


lim ^ 

13^00 V 


0 T ^ 

1 L = L>^. 


^See e.g., [O, §1.8]. A previous draft of the paper had stronger assumptions on 
E. We are grateful to a referee for pointing out to us that non-degeneracy suffices. 
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In particular, p^{x) is asymptotic in fl to the rational {d — 1)-cycle 
Consequently, in the low temperature limit^ the higher Boltz¬ 
mann distribution rationally quantizes. 

1.4. Improved higher matrix-tree theorems. Let X be a finite 
connected CW complex equipped with scalar structure i?*. Fix d < 
dimX. As above, we set E = Ed-i and W = Ed. 

Definition 1.18. Let \ C'd_i(X;R) —)■ be the linear 

transformation given by 

e-P^d*e^E^ 

where 

• d* : C'(i_i(X; R) —>■ (^^(X; R) is the formal adjoint to the bound¬ 
ary operator with respect to the standard inner product; 

• : Cd{X- R) ^ Cd{X- R) is given by 6 ^ for h e X^; 

• e^E. c^_^(X;R) ^ Cd_i(X;R) is given by / ^ e^^Ef j ^ 
Xd-i. 

The (restricted) biased Laplacian is the operator 

CE^^ ;= 5d-i(X;R) ^ 5d-i(X;R). 

Remark 1.19. The operator is just the formal adjoint to the 

boundary operator d with respect to the modified inner products {—,—)e 
and (—, —)w- 

The case E = was considered in [CCK] (in this instance we simplify 
notation and write for REW-^ similarly, if IF = 0 we write £^). 

Recall from [CCK, defn. 1.2] that a subcomplex T C X is a spanning 
tree (in dimension d < dimX) if 

. Hd{T) = 0 , 

• l3d-i{T) = j3d-i{X), and 

• C T C 

Remark 1.20. The above definition of spanning tree is a slight gen¬ 
eralization of the notion given in [CCK] in that we do not assume 
d = dimX. The definition given here is a matter of convenience only 
and doesn’t change any results of that paper since the spanning trees 
in dimension d of X are precisely the spanning trees of When 

dimX = d = 1, the definition is equivalent to the usual notion of 
spanning tree of a graph. 

For a finite complex Y and a choice of degree d < dim Y, let 
(7) dyGN, 

^The limit j3 ^ oo \s called the low temperature limit (cf. [CCK]). 
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be the order of the torsion subgroup of Hd-i(Y-,Z). Dehne the weight 
of a spanning tree T to be 

(8) Wt = Wt{W) := n • 

b£Ta 

Theorem C (Improved Higher Weighted Matrix-Tree Theorem). For 
a finite connected CW complex X and a degree d < dimX, we have 

L T 

where the first sum is indexed over all (d — l)-dimensional spanning co¬ 
trees, the second sum is indexed over all d-dimensional spanning trees, 
with Wt is as in (8) and tl as in (6). 

Remark 1.21. Theorem C effectively places spanning co-trees on the 
same footing as spanning trees. The special case E = t) recovers [CCK, 
thm. C], where the number fix appearing there is identihed here with 

It is worth singling out the special case E = = W\ 

Corollary D (Improved Higher Matrix-Tree Theorem). 

det C = ^ «i) ^r) ) 

L T 

where C = dd*: Hd-i(d^;M) —)■ ai is as in (5) and 9 t is 

as in (7). 

Remark 1.22. If the Betti numbers fij{X) are trivial for j = d — l,d — 2, 
then Corollary D reduces to the main result of [DKMl]. The reduction 
follows directly from the identity 

_ grf-i(X)grf_2(L) 

9d-2{X) ’ 

where 9j{Y) denotes the order of the torsion subgroup of Hj(Y;Z) 
(note: 9d-i{Y) is 9y appearing in (7)). The identity is an easy conse¬ 
quence of the short exact sequence of hnite abelian groups 

0 —)■ Hd-i{X) —)■ H(i-i{X, L) —)■ H(1-2{L) —)■ Hd-2{X) —)■ 0 , 

where homology is taken with integer coefficents. Here, we have used 
the fact that Hd-i{L) = 0, since fid-fiL) = fid-fiX) = 0 and dimL < 
d — 1. 

Outline. In §2 we develop the elementary properties of spanning co¬ 
trees. §3 contains the proof of the Boltzmann Splitting Formula (The¬ 
orem A). In §4 we introduce “tree-co-tree” duality which describes 
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a bijection between the spanning trees in a chain complex with the 
spanning co-trees in the dual cochain complex. In §5 we prove the 
Improved Higher Weighted Matrix-Tree Theorem (Theorem C). In the 
appendix (§ 6 ) we prove a summation formula for the pseudo-inverse to 
the boundary operator which unihes both the Kirchhoff projection and 
Boltzmann splitting formulas. 

Acknowledgements. We are indebted to a referee for valuable sugges¬ 
tions leading to a vastly improved version of the paper. Much of our 
writing was done while the last author was visiting the University of 
Copenhagen. He is indebted to Lars Hesselholt for providing him with 
support from the Niels Bohr Professorship. 

This material is based upon work supported by the National Science 
Foundation Grant CHE-1111350 and the Simons Foundation Collabo¬ 
ration Grant 317496. 


2 . Spanning co-trees 

If F is a held, recall that a fc-chain c G C'fc(W;F) is any F-linear 
combination of fc-cells. If 6 G Xk is a fc-cell, we write (c, h) for the 
coefficient of b appearing in c. If (c, b) 7 ^ 0, we say that b appears in c. 

Definition 2.1. A fc-cell b G is said to be essential if there exists 
a /c-cycle 2 ; G Q) such that {z^ b) 7 ^ 0. 

Lemma 2.2 ([CCK, lem. 2.2]). Adding or removing an essential d-cell 
from X increases or decreases (dd{X) by one, respectively, and fixes 
fid-iiX). 

Lemma 2.3. X has a spanning co-tree. 

Proof. The homomorphism Q) — Hd-fiX^Q) is surjec¬ 
tive with kernel Ki := Bd-i{X]Q). Set := Suppose that 

c G Bd-i{X] Q) is nontrivial. Let 6 be a (d — l)-cell of X that appears 
in c. Let be the result of removing b from The homomor¬ 

phism Hd-fiY^^^Q) Hd-i{X-,Q) is surjective; let be its kernel. 
Then the rank of is strictly less than that of by Lemma 2.2. 
Furthermore, fid- 2 iX‘^) = fid- 2 {X). By iterating (with Y^ replacing 
Y^, etc.) we eventually obtain a subcomplex Y^ C such that 

Hd-i{Y^] Q) —)■ Hd-i{X]Q) is an isomorphism. Then Y^ is a spanning 
co-tree. □ 

Proposition 2.4. Let ¥ be a field of characteristic zero. Let L (Z X 
be a {d — 1)-dimensional subcomplex that contains Then L is a 

spanning co-tree if and only if the composition 

(9) C',_i(L;F) ^ Cd_i(X;F) ^ C,_i(X; F)/H,_i(X; F) 
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is an isomorphism. 

Proof. Clearly, it suffices to prove the assertion when F = Q. Suppose L 
is such that (9) is an isomorphism. Consider the following commutative 
diagram: 

Z,_i(L;Q) ^ Z,_i(X;Q)-^. hfd-i(X;Q) 

k j 

Cd-imQ)^Cd-i{X-,Q)^Cd-i{X-,Q)/Bd-i{X-,Q). 

The left square is a pullback and the right square is a pushout. By as¬ 
sumption, the bottom composite is an isomorphism, so the top compos¬ 
ite is also an isomorphism. Therefore, iL*- Hd-i{L;Q) —)■ Hd-i{X]Q) 
is an isomorphism. The remaining two conditions of Dehnition 1.9 are 
easily verihed. Consequently, L is a spanning co-tree. 

For the converse, let x G (7^-1 (L;Q) be such that (tt o a){x) = 0. 
Then a{x) G Bd-i{X]Q) C Zd-i{X]Q). Since the left square is a 
pullback, we infer that x G Zd-i{L; Q). But p o is an isomorphism, 
and j is injective, so a: = 0. This establishes the injectivity of (9). 

For surjectivity, let G Q)/i?(i_i(X; Q). Lift this to any 

element y G Cd-i{X;Q). Then d{y) G (7^-2(T;Q) = (7^-2(X;Q) 
lies in Zd_ 2 (T;Q) since = 0. Furthermore, the pushforward of 
the homology class [c?(i/)] G Hd-2{L;Q) in Hd-2{X;Q) is trivial, since 
Hd- 2 mQ,) = Hd- 2 {X]Q). It follows that d{y) lies in Bd- 2 {X-,Q) = 
Bd- 2 {L-,Q). Hence, d{y) = d{x) for some x G Crf_i(L;Q). Then 
a{x) —y lies in Zd-i{X] Q), and since L is a spanning co-tree, there ex¬ 
ists x' G Zd-i{L] Q) so that 7 r(a(a;) —y) = (j opoi]f)(^x'). But z = vr(?/), 
so 

2 ; = Ti{y) = 7r(a(a;)) - i{p{iL{,x'))) = 7r(a(a; - k{x'))). 

We conclude that (9) is surjective. □ 

Remark 2.5. A referee has pointed out to us that Lemma 2.3 as well 
as Proposition 2.4 admit alternative proofs using matroids. 

Lemma 2.6. Let ¥ he a field. Then a splitting of the quotient homo¬ 
morphism Cd-i{X;¥) —)■ Cd-i{X-,¥)/Bd-i{X;¥) induces by restriction 
a splitting of the quotient homomorphism Zd-\{X]¥) Hd-i{X;¥). 

Proof. Consider the following commutative diagram, with exact rows. 

0-^ Bd-i{X- F)-. Zd-i{X- F)-^. Hd-i{X- F)-. 0 


0 


^Bd-i{X-¥) -.Crf_i(X;F) 


TT 


>C'd-i(X;F)/H,_i(X;F)->0. 





















14 


M. J. CATANZARO, V. Y. CHERNYAK, AND J. R. KLEIN 


Since Hd-i{X;¥) C Cd-i{X]¥)j we can restrict the given 
splitting to get a map Hd-i{X]¥) C(i_i(X;F). A simple diagram 

chase shows that this map factors through Zd-i{X]¥). □ 

Notation 2.7. For i < j and Y C X be a subcomplex. Set 

. 

Then j C j is a subcomplex of dimension < j. 

Corollary 2.8. Assume d > 2. Then the operation L ha Ld^d -2 defines 
a bijection between the spanning co-trees of X and the spanning co-trees 
of Xd 4-2 ■ Furthermore, = ai ■ 

Proof. By dehnition Cd-i{Xd 4 - 2 ]¥) = Zd-i{X]¥), so the diagram 
Cd-i{L-¥) ^ Cd-i{X-¥) ^ Cd-i{X-¥)/Bd-i{X-¥) 
coincides with the diagram 

Zd-i{Ld4-2]F) —^ Zd-i{Xd4-2]F) —^ Hd-i{Xd4-2]F) ■ 

It follows that the hrst part amounts to a restatement of Proposition 

2.4. 

To prove the second part, use the homotopy pushout diagram 

L -^X 

Ld,d-2 -^ Xd4-2 

and the long exact sequences in homology associated with the hori¬ 
zontal maps. Using the hve-lemma we infer that the homomorphism 
HfiX,L) —)■ HfiXd,d- 2 ,Ld,d- 2 ) is an isomorphism in all degrees. □ 

Remark 2.9. Corollary 2.8 reduces the proof of Theorem A to the spe¬ 
cial case of CW complexes X having trivial {d — 2)-skeleton (if d < 1, 
triviality is automatic). For such complexes the number ai coincides 
with the order of cokernel of the homomorphism 

Cd-imZ) ^ Cd-i{X;Z)/Bd-i{X;Z), 

The triviality of the {d — 2)-skeleton implies that the displayed map 
coincides with the homomorphism Hd-i{L]Z) —)■ Hd-i{X-,Z). 

Let Hd-i{X;Z)tor C Hd-i{X;Z) be the torsion subgroup. Let bi 
denote the order of the cokernel of the composite map 

(10) Hd-fiL-Z)^Hd-fiX-Z)^Hd-fiX-Z)/Hd-fiX-Z\,,. 
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Then (10) is a monomorphism of hnitely generated free abelian groups 
with hnite cokernel. Up to sign, the determinant of (10) is well-dehned 
and coincides with hi (see [CCK, prop. 6.63]). Furthermore, 

( 11 ) o,L = dxbL , 

where 6x is the order of ifrf_i(X; Z)tor- 

To complete the proof of Theorem A, we will construct a splitting 
of the map (7^-1 (X;M) —>■ Crf_i(X; M)/i?(i_i(X; R) that will give the 
relevant summation formula. For this, we shall appeal to the theory of 
generalized inverses. 

3. The proof of Theorem A 

3.1. Generalized Inverses. The theory of generalized inverses was 
developed to study linear systems Ax = h for which A~^ does not 
exist. Let A be an m x n matrix over R, and let b G R™ be given. 
Consider the linear system Ax = b. In general, such systems need not 
have a (unique) solution. One way to study the system is to attempt 
to minimize the norm of the residual vector Ax — b. Among all such x 
for which the norm of Ax — b is minimizing, we impose the additional 
constraint that the norm of x is minimizing. This is called a least 
squares problem.^ 

Remark 3.1. When A is surjective the residual vector having minimum 
norm is the zero vector. In this case the least squares problem reduces 
to the problem of hnding a solution of Ax = b such that the norm of x 
is minimized. 

The Moore-Penrose pseudo-inverse gives a preferred solution to 
the least squares problem. If b lies in the image of A, then a solution to 
Ax = b exists and the Moore-Penrose solution A'^b will be a solution 
having the smallest norm. Furthermore, the matrix A"*" exists and is 
unique [P], [BG, p. 109]. 

The operation A i—)■ A"*" satishes the identities 

(12) A+ = A*(AA*)+ = (A*A)+A*, 

where A* is the transpose of A (cf. [BG, chap. 1.6, ex. 18(d)]). In 
particular, when A is surjective, we obtain the formula 

(13) A+ = A*(AA*)-C 

^This is slightly more general than the usual formulation. The classical least 
squares problem assumes that A is injective. We will be primarily concerned here 
with the case when A is surjective. 
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Remark 3.2. If A is surjective, then one may drop the requirement that 
the target of A is based. That is, suppose more generally that A: M” —)■ 
1/ is a surjective linear transformation where V is not necessarily based. 
Then the least squares problem as well as the formula (13) make sense 
if we use the formal adjoint A*: V* —)■ (M"')* = in place of the 
transpose. Similarly, if A is injective, we may drop the requirement 
that the source of A is based. 


We will need a weighted version of the least squares problem. For 
this, we weight the standard basis elements of R"' by means 

of a positive functional p: —)■ R+. Then /i dehnes a modihed 

inner product (—, —)^ on R"', determined by {ei,ej)^ := ^{ei)Sij. The 
weighted least squares problem is to minimize \Ax — b\ such that \x\fj_ 
is also minimized. Again, the solution x = A'^h exists and is unique, 
where now A"*" is the weighted version of the Moore-Penrose pseudo¬ 
inverse. 

Assume now that A has rank m, i.e., A is surjective. For a subset 
S'c{l,...,n}of cardinality m, let be the mxm submatrix whose 
rows correspond to indices in S: 


{As)ij := Aij , for i = 1 ,... m, j e 5 . 

We will consider only those S such that A^ is invertible. Let is'- R”* —t 
R"’ denote the inclusion given by the rows corresponding to S. Set 

:= det(A5)2 JJ ^ 




F(e*) ’ 


and set V := state the summation formula for A+ 

in the case of surjective A. 


Theorem 3.3 (cf. [Be, Thm 1], [BT, th. 2.1]). Let A be an m x n 
matrix of rank m defined over R. Then the weighted Moore-Penrose 
pseudo-inverse of A is given by 

^ X]^5*5(As)“^ , 

s 

where the sum is taken over all indices S C {1, 2,..., n} such that As 
is invertible. 


Remark 3.4. Theorem 2.1 of [BT] concerns the case when A has rank 
n, i.e., when A is injective with arbitrary weights. The main result of 
[Be] applies to general A in the unweighted case /r = 1. When A is 
surjective, it is straightforward to deduce the weighted case from the 
unweighted one by the following transformation: replace A by A = 
AM“^, where M is the diagonal matrix having entries pi{ei) and 
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replace xhy x = Mx. This converts the weighted least squares problem 
{Ax = 6, p) to an equivalent unweighted problem {Ax = 6 , 1 ). The 
formula displayed in Theorem 3.3 is easily deduced from this, and we 
will omit the argument. 

Remark 3.5. Suppose that A : R” —>■ 1/ is a surjective linear transforma¬ 
tion in which V is not necessarily based (cf. Remark 3.2). Furthermore, 
suppose that FT C R is a lattice, i.e., a hnitely generated abelian sub¬ 
group such that the induced map FT 0z R —)■ R is an isomorphism. 
Then a choice of basis for H determines one for R and Theorem 3.3 
applies. 

Moreover, the formula is invariant with respect to base changes for 
Ff, since the numbers ts are squares of determinants. Consequently, 
Theorem 3.3 is really a statement about surjective linear transforma¬ 
tions A: —)■ R for vector spaces R that come equipped with a 

preferred lattice H. 

Proof of Theorem A. By Remark 2.9 there is no loss in assuming that 
X has trivial (d—2)-skeleton. By Remark 1.8 and Lemma 2.6, it suffices 
to produce a splitting of the quotient homomorphism vr: Cd-i{X] R) —)■ 
Cd-i{X ; R) / Bci-i{X ; R). Here we use the weighted basis of Cd-i{X ; R) 
defined by the cells and the weighting given by 6 i—)■ Use the lattice 

H C Cd-i(X; R)/i?d-i(df ; R) given by the image of the homomorphism 
Cd-i{X-,Z)/Bd-i{X-,Z) ^ C',_i(X;R)/R,_i(X;R). 

Applying Theorem 3.3 and Remark 3.5 to tt gives a splitting, writ¬ 
ten as a sum over subsets S of the basis elements of Cd-i(X;R). By 
Proposition 2.4, the collection of these subsets are in bijection with the 
set of spanning co-trees. The inclusion is corresponds to the inclusion 
Crf_i(F;R) —)■ C'd_i(X;R) and corresponds to As- Then the deter¬ 
minant of As is given up to sign by := tl/9x (cf. Remark 2.9) and 
the prefactor is given by the reciprocal of V' := Hence, the 

desired splitting is given by 

Remark 3.6. If we fix a function W : Xd —)■ R, we may instead apply 
[Be, Thm 1] to the inclusion map g: Zd(X;R) —>■ C'd(X;R). This 
produces an orthogonal splitting Crf(X;R) Zd{X-,'R) to q in the 
modified inner product on (^^(X; R). The splitting is written as a sum 
indexed over the set of spanning trees as in [CCK]. In fact, this gives 
quick alternative proofs to Theorem A and Addendum B in [CCK]. 
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4. Tree-co-tree duality 

Let X be a finite connected CW-complex and fix a natural number 
d < dimX. We let 

{C,{X),d) and 

denote the cellular chain and cochain complexes of X over Z. The set 
Xd determines preferred isomorphism C"^(X) = Cd{X). With this iden¬ 
tification, the coboundary operator S: C‘^~^{X) C"^(X) corresponds 

to d*: Cd-i{X) —>■ Cd{X), the formal adjoint to d (here, coefficients 
can be taken in any commutative ring). 

Denote by Xd{X) and XJ(X) the sets of d-dimensional spanning 
trees and spanning co-trees, respectively, of X. For i < j, recall the 
subquotient Xj^i = X^-^^X^. 

Lemma 4.1. There are preferred hijections 

Xd{X) = Xd{Xd,d- 2 ) and X^i(X) = X*d_,{Xd,d- 2 ) • 

Proof. The second bijection is merely a restatement the first part of 
Corollary 2.8. For T G Xd(X) and L G Jy[_i(X), the assignments 

T HA Td,d-2 and L ha Ld,d-2 

define maps Xd(X) -)■ Xd{Xd,d- 2 ) and X2-i{X) -)■ Xl_i{Xd,d- 2 )- In¬ 
verse maps are defined as follows: given a spanning tree T' for Xd^d- 2 , 
we form a complex T by attaching the set of d-cells appearing in T' 
to It is straightforward to verify that T is a d-dimensional 

spanning tree for X. 

Similarly, if L' is a spanning co-tree for Xd^d- 2 , then the complex L 
given by attaching the (d —l)-cells of L' to gives a d-dimensional 

spanning co-tree. □ 

For T G Xd(X), let 6t be the order of the torsion subgroup of 
Hd-i{T]Z) (cf. [CCK, p. 3]) and for L G XJ_]^(X), recall that ql 
is the order of Hd-i{X, L; Z). 

Lemma 4.2. ForT G Xd{X) and L G Jy[_i(X) we have 

and ai = ■ 

Proof. The first part follows immediately from the exactness of the 
sequence 

0 ^ Hd-i{T-Z) ^ Hd-i{Td,d-2]I^) ^ Lfd-2(T(''-2);Z) 

and the fact that Z) is free abelian. The second part is 

just a restatement of the second part of Corollary 2.8. □ 
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Remark 4.3. In view of these lemmas, all the relevant properties of d- 
dimensional spanning trees and {d — l)-dimensional spanning co-trees 
depend only on the boundary operator d : Cd(X) —)■ Cd-i(X) and are 
therefore reduced to properties of two-stage chain complexes of hnitely 
generated free abelian groups, i.e., to statements in linear algebra over 
Z. We will now make this precise. 

For a hnite set S, let be the free abelian group with basis set S. 

Definition 4.4. For hnite sets P and Q, set A = Z^ and B = Z^. 
Let d \ A ^ B he a. homomorphism. A spanning tree of d consists of a 
subset T G P such that the composition 

Z^ ^ A^B 

is an isomorphism modulo torsion. The set of spanning trees of d is 
denoted hy P{d). 

Similarly, a spanning co-tree of ci is a subset L G Q such that the 
composition 

(14) Z^ gB^ B/d{A) 

is an isomorphism modulo torsion. The set of spanning co-trees of d is 
denoted by iF*{d). 

For an abelian group U let U* := homz{U,Z). This dehnes an con- 
travariant endo-functor on abelian groups. Let d* \ B* — )• A* be the 
homomorphism induced by d. For a subset S' C P, let S'-*- = P \ S' 
denote its complement. The proof of the following is a straightforward 
exercise left to the reader. 

Lemma 4.5 (Tree-Co-tree Duality). The operation T ha T-^ induces 
a bijection P{d) = iF*{d*). 

4.1. Finite chain complexes. A Z-graded chain complex C* over Z 
is finite if it is degree-wise hnitely generated and free and has hnitely 
many non-zero terms. If C* is hnite then so is its linear dual 

PC'* = hom(C_*, Z). 

If we hx d G Z, then the notion of d-dimensional spanning tree and co¬ 
tree is dehned in this context. Let Pd(C*) be the set of d-dimensional 
spanning trees of C*. Similarly, we let P))(C*) be the set of d-dimensional 
spanning co-trees of C*. The following is an immediate consequence of 
Lemma 4.5 combined with Remark 4.3. 

Corollary 4.6 (Chain Tree-Co-tree Duality). Assume C* is finite. 
Then for d G Z, there is a preferred bijection Pd(C*) — XffiDCfij. 
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Remark 4.7. It is not difficult to give a spectrum-level version of Corol¬ 
lary 4.6. When X is a finite CW spectrum, the spectrum DX (corre¬ 
sponding to the linear dual of a chain complex) is the Spanier-Whitehead 
dual of X, i.e, the function spectrum F(X, S'°). 

Corollary 4.6 may be related to some of the results of [MMRW]. 

5. The proof of Theorem C 

Given a scalar structure on X and a dimension d < dimX, one has 
a pair of biased Laplacians defined by the commutative diagrams 
(15) 

Bd-i (X; R) B\X- R) B\X] R) 



Bd-i{X- R) R'^(X; R) R). 


Observe that 

• 

The operators invertible and have the same de¬ 

terminant. To avoid notational clutter, when E and W are understood 
we simplify notation and set C := and C* := Cw,e- 

Recall that the goal is to exhibit a decomposition of the determinant 
of T as a sum over trees and co-trees: 

(16) det(£) = det(£*) = ^ ^ t'lWt ■ 

We start with a weak version that establishes (16) up to a factor 
that does not depend on either E or W. We consider T as a function 
of the pair (T, W) taking values in the space of linear operators on 
Rrf_i(X; R). We then measure the variation of In det(T) in W assuming 
E is held hxed. This gives 



dlndet(/:) = -Tr(ddWe-^d*e^£~^) 

(17) 

= -Tr(dWe-^d*e^£-^d) 


= -Tr(dWA ), 

where 


(18) 

A := Bd-i(X;R) Cd(X;R) 


(compare [CCK, eqns. (6),(7)]). It is straightforward to check that A 
is the pseudo-inverse for d : Cd(X) —)■ Bd-i(X). Thus we can follow. 
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mutatis mudandis, the proof of [CCK, prop. 4.2] (which is the E = 0 
particular case of C) to arrive at 

(19) det(£) = det(£*) = j{E) wt{W) . 

T 

where the prefactor ^{E) is independent of W. 

We next make use of duality, i.e., det(£) = det(£*) = d*e^de~^, 
so that variation over E, with W held hxed, is equivalent to the case 
considered above. Consequently, 

( 20 ) det(/:) = det(£*) = TiW) tl{E) , 

L 

where the prefactor 7 *(hh) is independent of E. 

Setting (19) and (20) equal, we immediately obtain the following 
weak form of Theorem C: 

(21) det(£) = det(£*) = 7(^rL)(^M;T), 

L T 

in which the prefactor 7 is a constant independent of E and W. The 
proof of Theorem C is complete once we establish the following claim. 


Claim. 7 = 1/6'!' 


To prove the claim, by the argument of [CCK, § 6 ] it is enough to 
restrict to the case when X is a spanning tree of dimension d. Fur¬ 
thermore, by Lemmas 4.1 and 4.2 we can assume that X = Xd^d -2 = 
XC)/xC- 2 )_ gince 7 is independent of E and W we can further assume 
that E = 0 = W. In this instance, by [CCK, cor. D] it will suffice to 
establish the identity 

(22) /ix = ^ ai) 

L&T*(X) 


where fix denotes the square covolume of the lattice Bd-i{X]'L) C 

Bd-i{X;R). 

With respect to our hypotheses, the cellular chain complex of X 
is determined by the (rationally injective) homomorphism d: A ^ B 
of hnitely generated free abelian groups. By slight abuse of notation, 
denote this chain complex by X and let Y be the dual chain complex 
given by d*: B* ^ A*, where A* = hom(y4, Z). Then 


(23) 


det(£) = det(£*) 
= det(c}*9) 


_ t^Y 






T ■ 


T&T{Y) 
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where the last equality follows from [CCK, cor. D], In the above, we 
have implicitly identihed d* with the transpose of d using the preferred 
bases. The number /iy is the square of the covolume of the lattice 
dehned by the image d*{B*) C d*{B*) (g)^ R where the latter term is 
given an inner product by declaring it to be an isometric subspace of 
A* M = hom(74, R). By Lemma 4.5, B{Y') = B*(X). 

Given a spanning co-tree L G let denote the order of the 

cokernel of L —)■ B/d{A). For a subgroup T C A let 6t denote the order 
of the torsion subgroup of the cokernel of the composition T ^ A ^ B. 
Note that these definitions coincide with the ones we gave in the case 
of CW complexes with trivial {d — 2)-skeleton. 

The identity (22) (and hence the claim) is now an immediate conse¬ 
quence of the following: 

Lemma 5.1. With respect to the above assumptions, we have 

(1) det(£) = /ix; 

(2) 6y = Ox; 

(3) /iy = Oy; 

(4) ol = 

Proof. Statement (1) is a special case of [CCK, cor. D]). 

Let C be the cokernel d. Then we have a short exact sequence 
0—>-(7—)-0. Applying hom(—,Z) gives a short exact 
sequence 

^ B* ^ A* ^ ext(C, Z) ^ 0. 

Since C is hnitely generated, the group ext(C', Z) is a torsion group 
which is (non-canonically) isomorphic to the torsion subgroup of C. 
This gives (2). 

Statement (3) follows from the discussion after dehnitions 6.2 and 
6.5 of [CCK] as applied to the real isomorphism d*{B*) —)■ A*. Lastly, 
(4) is easily deduced from the fact that the cokernels of Z^ —)■ 5/ d{A) 
and A —>■ Z^^ are canonically isomorphic. □ 

This establishes the claim and completes the proof of Theorem C. 

Remark 5.2. There is different proof of the claim that makes use of the 
“low-temperature limit” method of [CCK, §5] (with respect to E and 
W) as applied to equation (21), with the sum on the right-hand side 
dominated by one term. However, the proof we have given here has 
the advantage of being shorter and less technical. 
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6. Appendix: Kirchhoff-Boltzmann unification 

For any fixed dimension d, the bonndary operator d: C'rf(X;M) —)■ 
Crf_i(A;R) factors as 

in which the hrst map is surjection and the second is an injection. 

With respect to these inner product structures, one infers that the 
pseudo-inverse of d: —)■ C'(i_i(X;R) is given by the composi¬ 

tion 

C'd-i(A;R) ^ B,_i(A;R) ^ Cd{X-R) 
consisting of the pseudo-inverses of i and p. (cf. [BG, p. 48, ex. 17]). 
Here, is the pseudo-inverse of p with respect to the modified inner 
product structure on C'd(A;R) (cf. Remarks 3.2 and 3.5). Similarly, 
is the pseudo-inverse for i defined using the modified inner product on 
Grf_i(A;R). 

Therefore, all three pseudo-inverse formulas will follow from one, e.g., 
from the formula for the surjection, since the formula for an injection 
can be obtained by taking transposes and applying duality (Lemma 
4.5), whereas the general pseudo-inverse is obtained by taking the com¬ 
position. 

For a tree T G J7i(A), the composition 

Cd{T-,M.) ^ C'd(A;R) 4 Bd-i(A;R) 

is an isomorphism. Let ax denote its inverse. Let px '■ -Sd-i(A;R) 
Grf(A;R) be the composition 

Rd-i(A;R) ^ C,{T-R) 4 C^^X-R ). 

Then using Theorem 3.3 one has 

(24) Vw = Aw-.= '^wx, 

T T 

where the sum is over spanning trees T G Xd{X). 

Similarly, the pseudo-inverse of the inclusion i : Bd-i{X]M.) —)■ 
Grf_i(A;R) is obtained as follows: for a co-tree L G X^_i{X) the com¬ 
position 

B,_i(A;R) ^ G,_i(A;R) ^ G,_i(A;R)/C',_i(L;R) 

is an isomorphism; let (3x denote its inverse. Let (x ■ 4-1 (A; R) 
4_i(A;R) be the composition 

4-i(A;R) ^ 4-i(A;R)/4-i(i:;K) A 4_i(A;R) 
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where the first map is vector space projection. 

Then 

(25) 4 = VE = ^ri, 

L L 

where the sum is over spanning co-trees L G 

Combining (24) and (25) with the pseudo-inverse composition prop¬ 
erty we arrive at a formula which encompasses both the Boltzmann 
splitting formula (Theorem A) and the higher Kirchhoff projection for¬ 
mula of [CCK, thm. A]: 

Theorem 6.1 (Kirchhoff-Boltzmann Projection Formula). The pseudo¬ 
inverse of the boundary operator d: Crf(X;M) —)■ Cd_i(X;M) with re¬ 
spect to the modified inner products defined by E and W is given by 

(26) ^LWTaL,T , 

L,T 

where aL,T = Tt°Cl ■ C'(i_i(X;M) —>■ C(i(X;M) and the sum is indexed 
over L G iF^_fiX) and T G iFd{X). 

Remark 6.2. Theorem 6.1 gives a concrete expression for the average 
current in the case of periodic stochastic driving (in the adiabatic limit) 
[CKS]. Let 7 be a smooth 1-dimensional cycle in the vector space of 
parameters {E,W) and let x G Zrf_i(X;Z) be a (d — l)-cycle. Let 
X = [x] G ifrf_i(X;Z) be the associated homology class. Then the 
average current of ( 7 , [x]) is dehned to be the d-cycle 

q-.= f Adp^ e Zd{X-,R), 

J-f 

where A = A(E, W) = is the operator of equation (18) 

and p® is the higher Boltzmann distribution of x (see [CKS] if d = 
dimX = 1 and [C] for d > 1.). Then, after some straightforward 
algebraic manipulation using Theorem 6.1, we hnd 

q = / dgfi A dpl , 

L,T 

where gf = gfi(W) := WT/Xx^gf = gf^E) := Tl/Xl and a is any 
smooth 2 -dimensional chain in the space of parameters satisfying da = 

7- 
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